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TECHNICAL NOTE 2140 

THEORETICAL ANTISYMMETRIC SPAN LOADING FOR WINGS OF 
ARBITRARY PLAN FORM AT SUBSONIC SPEEDS 
By John DeYoung 


SUMMARY 


A simplified lifting-surface theory that includes effects of com- 
pressibility and spanwise variation of section lift-curve slope is used 
to provide charts with which antisymmetric loading due to arbitrary anti- 
symmetric angle of attack can be found for wings having symmetric plan 
forms with a constant spanwise sweep angle, of the quarter-chord line. 
Consideration is given to the flexible wing in roll. .Aerodynamic char- 
acteristics due to rolling, deflected ailerons, and sideslip of wings 
with dihedral are considered. Solutions are presented for straight- 
tapered wings for a range of swept plan forms. 


INTRODUCTION 


Reference 1 has been for many years the standard reference for esti- 
mating the stability and control characteristics of wings. The lifting- 
line theory on which this work was based gave generally satisfactory 
results for straight wings having the aspect ratios considered; however, 
the use of wing sweep combined with low aspect ratio has made an exten- 
sion of this work desirable. Lifting— line theory cannot adequately 
account for the increased induction effects due to sweep nwri low aspect 
ratio; consequently, it has been found necessary to turn to the more 
complex lifting— surf ace theories. 

Of the many possible procedures, a simplified lifting-surface theory 
proposed by Weissinger and further developed and extended in reference 2 
has been found especially suited to the rapid computation of character- 
istics of wings of arbitrary plan form. Comparisons with experiment have 
generally verified the theoretical predictions. In references 3 and 4, 
this method has been used to compute for plain, unflapped wings, the aero- 
dynamic characteristics dependent on symmetric loading. The same simpli- 
fied lifting-surface theory can be extended to predict the span loading 
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resulting from antisymmetric 1 distribution of the wing angle of attack. 
From such loadings, the damping moment due to rolling, the rolling 
moment due to deflected ailerons, and the rolling moment due to dihedral 
angle with the wing in sideslip can be determined. A recent publication 
(reference 5) makes use of the simplified lifting— surface theory to find 
span— loading characteristics of straight— tapered swept, wings in roll and 
loading due to dihedral angle with the wing in sideslip. Experimental 
checks of the theory for the damping— in— roll coefficient and rolling 
moment due to sideslip were very favorable. The range of plan forms 
considered in reference 5 is somewhat limited and aileron effectiveness 
was not included. The loading due to aileron deflection normally involves 
excessive labor when computed by means of the simplified lifting— surface 
theory; however, development of the theory, presented in reference 6, 
that deals with flap and aileron effectiveness for low-aspect— rat io wings 
provides a means by which the simplified lifting— surface method can be 
used to obtain spanwise loading due to aileron deflection. 

It is the purpose of the present analysis to provide simple methods 
of finding antisymmetric loading and the associated aerodynamic coeffi- 
cients and derivatives for wings with symmetric plan forms limited only 
by a straight quarter— chord line over the semi span. Means will be pre- 
sented for finding quickly the aerodynamic coefficients of span loading 
due to rolling, of span loading due to deflected ailerons, and of span 
loading due to sideslip of wings with dihedral. Flexible wings, when 
the flexure depends principally on span loading as in loading due to 
rolling, can be included in the analysis. 


NOTATION 


aspect ratio ( — 
v S' 


wing span measured perpendicular to the plane of symmetry, feet 
2 wing chord, feet 
c a 2 aileron chord, feet ' 
c av 2 mean wing chord ^ , feet 


local lift coefficient 


0 


local lift 'N 


qc 


J 


iThe word antisymmetric is understood to indicate that a distribution of 
loading or angle of attack is equal in absolute magnitude on each half 
of the wing but of opposite sign. 

2 Measured parallel to the plane of s y mm e try. 
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PDi 

Cl 

Cl 

4 P 

c lb 

cic . 
C l c av 
d v 

e nk 

G 

G 

G 


Induced drag coefficient ^ induced drag ^ 


rolling-Hiioinent coefficient . 


^rolling moment ^ 


rolling moment due to rolling 


qSb 

Sc i 


S(pb/2V) 




/ 

per radian 


/ SC Z \ 

rolling moment due to aileron deflection ( gg— J , per radian 

spanwise loading coefficient for unit rolling moment ^ ~q^ ) 
scale factor 

factors of loading interpolation function 
spanwise loading coefficient or dimensionless circulation 

Gfi) « GO 

spanwise loading .coefficient due to rolling ( ■ ; 

spanwise loading coefficient due to aileron deflection ^"g- ^ } 
per radian 

wing geometry, compressibility, and section lift-curve— slope 
parameter [ d v(^)(^7f)_ 

integration factors for spanwise loading due to ailerons 
Mach number 


M 


per radian 


m 


arbitrary number of span stations defined by tj = cos 

m+1 


p rate of rolling, radians per second 

pb/2 V wing— tip helix angle, radians 

Pyu coefficient depending on wing geometry and indicating the 

influence of antisymmetric loading at span station n on 
the downwash angle at span station V 
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q free— stream dynamic pressure, pounds per square foot 

S wing area, square feet 
t 2 ratio of aileron chord to wing chord 
V free— stream velocity, feet per second' 

w induced velocity, normal to the lifting surface, positive for 
downwash, feet per second 

y lateral coordinate measured from the wing root perpendicular to 
the plane of symmetry, feet 

a v 2 section angle of attack at span station .v, radians 

£a v 2 angle of antisymmetric twist of the elastic wing produced by the 
loading due to rolling, radians 

da - 

— rate of change of wing— section angle of attack with control— sur— 
d-S face angle for constant section lift coefficient 


3 angle of sideslip, radians 

T dihedral angle measured perpendicular to the plane of symmetry, 
radians 

r c spanwise circulation, feet squared per second 

5 2 angle of deflection of full wing— chord control surface, radians 

8 angle of deflection of full-wing— chord control surface, measured 

perpendicular to the hinge line, radians 

q dimensionless lateral coordinate 

q a dimensionless aileron span 

^c.p. spaE-wise center of pressure on one wing panel 

0 trigonometric spanwise coordinate 9, indicating the edge of the 
aileron span, radians 

2 See footnote 2 p,2. 



aileron span 

b7£ 




compressibility parameter 
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K 


V 


A 

A 

9 


n,V 


k 

c.p. 

a 

t 

T 

E 

av 


ratio of section lift-curve slope at a span station v to — , 
both at the same Mach number ^ 

sweep angle of the wing quarter— chord line, positive for sweep- 
back, degrees 


compressibility sweep angle parameter 

taper ratio ( tip chord \ 

\ root chord J 


tan' 


_! / tan A 

V pa 


degrees 


trigonometric spanwise coordinate (cos^ 1 tj), radians 


Subscripts 


integers pertaining to specific 


T) = cos 


nn 

8 


V JT 

or tj = cos ^ 
8 


span stations given by 


pertaining to span station k 

center of pressure 

aileron 

pertaining to fraction-of -wing-chord ailerons 
wing tip 
wing root 
average or mean 


DEVELOPMENT OF METHOD 


The simplified lifting— surf ace method used herein replaces a lift- 
ing surface by a lifting vortex located at the wing one-quarter-chord 
line. The boundary condition for determining the vortex strength dis- 
tribution specifies that, along the three— quarter— chord line of the 
wing, there shall be no flow through the lifting surface. In 'effect, 
this specifies that, at the three— quarter— chord line, the ratio of the 
velocity normal to the mean camber line (induced by the bound and trail- 
ing vortices) to the velocity of the free stream shall equal the sine 
of the angle of attack. 
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Span loadings are theoretically additive. Since the symmetric 
angle-of-attack distribution contributes only to symmetric loading, it 
follows that the antisymmetric loading is independent of symmetrically 
distributed wing twist or camber; hence, to find antisymmetric loading, 
it is only necessary to consider the loading resulting from the anti- 
symmetric distribution of the angle of attack across the wing span. In 
the subject case, such a distribution is experienced by the wing as 
induced angle due to rolling , the effective twist due to aileron 
deflection, or sideslip of the wing with dihedral. 

For an antisymmetric angle-of-attack distribution, the loading dis- 
tribution will be equal in absolute magnitude oh each semispan, but of 
opposite sign. The loading therefore needs only to be found over the 
semispan, and, since the loading is zero at the wing root, only span 
stations outboard need be considered. The mathematical development of 
the simplified lifting— surface method for. the case of antisymmetric 
loading is given in Appendix A. As shown in Appendix A, (m-l)/2 linear 
equations in terms of loading distribution are obtained which satisfy 
the wing angle-of-attack conditions at the three-quarter-chord line at 
m stations n, where m is an arbitrary odd integer. These equations 
are represented by the summations 

m-1 

2 

a-v = ^ Pva 

n=l 

^n considering the case of the angle induced by rolling as equivalent 
to an antisymmetric distribution of twist, it must be noted that 
account should be taken of the fact that a rolling wing leaves a 
twisted vortex trail; whereas a twisted wing does not. The difference 
in induction effects on the wing of the straight and twisted vortex is 
considered insignificant here, as has been assumed in other analyses. 

The reader should note that the boundary condition is given by 
w v = V sin a v from which (w/v)v is seen to equal sin ay. The sub- 
stitution of a v for sin a v has the effect of increasing the value 
of loading on the wing above that necessary to satisfy the boundary 
condition. However, the boundary condition was fixed assuming that the 
shed vortices moved downstream in the extended chord plane. A more 
realistic picture is obtained if the vortices are assumed to move’ down- 
stream in a horizontal plane from the wing trailing edge. It can be 
seen readily that, if this occurs, the normal component of velocity 
induced by the trails at the three— quarter— chord line is reduced and, 
if the boundary condition is to continue to be satisfied; the strength 
of the bound vortex must increase. It follows that substitution of a v 
for sin cl\i then has the effect of accounting for the bending up of 
the trailing vortices. It is not known how exact the correction is, 
but the calculations and experimental verification show it to be of 
the correct order. 


v - 1, 2 , 3* • 


m-1 

2 


( 1 ) 
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where 

a V antisymmetric angle of attack at wing station v 

P Vn coefficients that for a given value of m depend on wing 
geometry, compressibility, and section lift-curve slope 

G n . loading coefficients at span stations n 


The application in Appendix A of the present report is with m=7 
Since the loading at the midspan station is known to be zero, considera- 

req ^ ired of on ^ V^ ee Nations: n=l,2,3, equal to wing semispan 
positions of q = cos ( nn/8) = 0.924; 0.707; and O. 383 . Equation (l) 
thus becomes ' ' 



n=l 


V = 1,2,3 


( 2 ) 


where the integer v pertains to span station tj 


cos (vjr/8) 


To obtain the loading coefficients 
only to evaluate the coefficients p^ 
the antisymmetric angle of attack a v . 


^n - ( c Z c /2t>) n , it remains 
and the spanwise variation of 


Evaluation of the Coefficients p Vn 


Since ^ m is chosen, p^ becomes a function only of wing geometry, 
compressibility, and section lift-curve slope. The effects of compressi- 
bility and section lift-curve slope are equivalent to a change in wing 
plan form and can be accounted for by. a proper adjustment of the p Vn 
values. As shown in Appendix B, p Vn can be conveniently presented as 
a function of two parameters, namely, a compressible-sweep-angle parame- 
ter defined as Ap= tan— !( tan A/p) and a parameter Hy involving the 

ratio of wing span to wing chord and variable section lift-curve slooe 
defined by * ’ 

H v ■ d v(^r) (jys) (3) 

where 


^ V ratio of experimental section lift — curve slope at span station v 
to the theoretical value of 2n/p, both at the same Mach number 


Compressibility and section lift— curve slope are discussed in Appendix B. 
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cy wing chord at span station V 

The value dy is a scale factor given by 

d v = 0.061 for v =1 
= 0.23^ for v =2 
"= 0.381 for v =3 


} W 


Equation (3) can be written in alternative form that gives Hy in 
terms of wing geometry parameters that are more significant; thus 


Hv ’ : K^)[ 


(Ky/K av ) ( C v / c av) 


(5) 


where 


av 

KV/^av 

c y/ c av 


ratio of average section lift-curve slope to 2«/p both at 
the same Mach number 

spanwise distribution of section lift-curve slope for a 
given Mach number 

spanwise distribution of the wing chord 


( PA/K av ) compressible aspect ratio and average section lift-curve— 
slope parameter 

The term —7 of equation (5) gives an effective aero- 

'Ky/K a -y) ( C y/ C av ) 

dynamic taper of a wing. The distribution of Ky/K av may vary with 
Mach number, particularly at transonic speeds (e.g., due to spanwise 
variation of airfoil section). However, since the distribution con- 
tributes to taper effect, the loading distribution and not the total 
loading will be appreciably affected. 


With Hy determined from equations (3) or (5) and (4), the values 
of Pvn> nine in all, are presented in figure 1 where pyjj is given 
as a function of Hy for various values of . 


For the case of straight— tapered wings with arbitrary section lift- 
curve— slope distribution for which the chord distribution is specified 
by taper ratio, evaluation of equation (5) is given in figure 2 where 

(tty/ Kay for each of the three span stations is shown as a function of 
( PA/K av ) 
taper .ratio. 
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Evaluation of Antisymmetric Angle-of-Attack Distribution a v 


The antisymmetric angle-of-attack distributions most commonly 
encountered are those resulting from rolling wings, aileron deflection, 
and sideslip of wings with dihedral. Evaluation of the angle-of-attack 
distributions for these various cases is outlined in the sections imme- 
diately following. 


Bolling wings .- For the case of the rigid wing, the induced veloc- 
ity nor ma l to the wing surface is equal to the upwash velocity experienced 
by the rolling wing. Thus, at span station V 



(6) 


where pb/2V is the tip helix angle. It should be noted that the rela- 
tion given by equation (6) assumes the wing structure to be rigid in that 
the distribution of a>y is completely defined by the linear distribution 
of helix angle. In the case of flexible wings, however, the expression 
for a y must be modified to account for the streamwise angle-of-attack 
change which may occur due to bending or torsional deflections. In this 
case, . ' 


a 


V 



+ Aa v 


(7) 


where Ax v represents the modifying influence of flexibility. Normally, 
is not considered for straight wings since only the effect of tor- 
• sion (which is usually small) is involved. On swept wings, however, the 
effect of bending can cause Ax v to be quite large so that the. a v 
distribution may be affected considerably. Due to the interaction 
existing between the aerodynamic and structural forces, Ax cannot be 
determined directly, but must be found through equations of ^equilibrium 
or by iteration. With the loading for the rigid wing provided, however, 
the iteration procedure becomes relatively easy to apply. The first 
approximation of a v is found from the loading of the rigid wing and 
further refinements of a v may be found utilizing the successive loadings 
for the flexible wing as determined. 


Deflected ailerons .— Where the spanwise distribution of the angle 
a v is to be considered equivalent to antisymmetric aileron deflection, 
it must suffer a discontinuity at the spanwise end of the control surface. 
The loading when such a discontinuity is present can be duplicated by a 
proper distribution of antisymmetric twist. In Appendix C, the anti- 
symmetric twist distribution required- by the present theory to give accu- 
rate span loading distribution due to ailerons is found with the aid of 
zero-aspect— ratio wing theory given by reference d. To minimize the 
confutation involved, it is convenient to consider both the case of out- 
board and inboard ailerons. 
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1. Outboard ailerons : With m=7, three different aileron spans 

can be conveniently defined for the outboard ailerons. For 
the aileron spans q a , measured from the wing tip inboard, 
the antisymmetric twist distribution required per unit deflec- 
tion of full-wing-chord ailerons, oy/6, is given by 


Case 

I 

II 

III 


0.169 

0.444 

0.805 

0-1 

b 

1.003 

0.971 

0.998 

Qg 

6 

.017 

• 996 

.991 

Oo 

6 

.006 

.014 

.978 


(8) 


2. Inboard ailerons : With m=7, three different aileron spans 

can be conveniently defined for the inboard ailerons. For 
the aileron spans q a measured from the wing midspan out- 
board, the antisymmetric twist distribution required per unit 
deflection of full-^wing-chord ailerons, ay/b, is given by 


Case 

IV 

V 

VI 

Ha 

0.556 

0.831 

1.000 

a-i 

b 

0.044 

0.013 

1.016 

9 

-.017 

.961 

.979 

n 

1.087 

, 

1.095 

1.101 


(9) 


Sideslip of wings with dihedral .— For calculating the rolling 
moment caused by dihedral angle for the sideslipping wing, the effect of 
the skewness of the vortex field in altering the effects of the dihedral 
angle will be assumed to be small (as assumed in reference 5)- The 
problem then simplifies to finding the rolling moment due to anti symmetric 
single of attack .with the unskewed vortex field. The solution to this 
problem is the same as for the. ailerons which has already been solved. 
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The antisymmetric distribution of angle of attack for the side- 
slipping wing with dihedral is gi T jn by 

4r (10) 

where 

cuy effective angle-of-attack distribution 

3 angle of sideslip measured positive in the counterclockwise 

direction from the plane of symmetry 

T* dihedral angle 

The wing parameter r ,i s not affected by compressibility. Equation (10) 
is approximate for small values of 3 and F. 

For unit 3 T over the span of the ailerons considered, 

P T = 5 (11) 

can be substituted for 5 in equations (8) and (9)* 


APPLICATION OF METHOD 


- For the cases of antisymmetric angle-of-attack distributions result- 
ing from rolling, aileron deflection, or sideslip with dihedral, it is 
possible to present a set of simultaneous equations which are required 
for the solution, of the load distribution for an arbitrary plan form. 

With the loading known, integration formulas can be given to determine 
aerodynamic coefficients. 

The loading-distribution coefficient Gn determined from the solu- 
tions of the simultaneous equations, are functions of p vn which has 
been shown in a preceding section to be a function of wing geometry, 
compressibility, and section lift-curve slope. The aerodynamic coeffi- 
cients are integrations of the load distribution and, therefore, will 
also be a function of wing geometry, compressibility, and section lift- 
curve— slope parameters. Application of the method to the general solu- 
tion for arbitrary chord distribution is outlined and solutions are pre- 
sented for the case of straight taper. 
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General Solution 


Aerodynamic characteristics due to rolling .— The solutions for the 
aerodynamic effects due to the rolling wing will be found and loading, 
rolling moment, spanwise center of pressure, and induced drag will be 
obtained . 

1. Simultaneous loading equations : The p^ values are 

obtained from figure 1 and table I s with values of Hv 
given by equations (3) or (5). 

The simultaneous equations (2), for the rigid and flexible 
wing, respectively, become: 


-0.924 = PnGx + Pi^* 2 + P 13 G 3 
-0.707 = P21 G 1 + P22 G 2 + P23 G 3 
—O.383 = PaiG-L + P 32 G 2 + P33 G 3 


where 


G n - 


pb/2V 


-0.924 + 


55727 = PllGl + ^ + P 13 G = 


-0.707 + ~ ^ 2lGl + ^22 G 2 + P23 G 3 


- AOq _ _ _ 

^• 3 3 + = P 31 G 1 + + P 33° : 


where G n = and i ncr e men 'tal angle of 
attack due to aeroelastic effects. 

Loading distribution : The loading-distribution coefficient 
is given by G = Cjc/2b. Other forms of the loading coeffi- 
cient are given by the identities 


3 Values of p^ beyond the scope of figure 1 are included in table I. 
For values of Hy larger than those included in figure 1 and table I, 
the p^ curves can be obtained from equation (B8) which gives the 
linear asymptotes of the p^ function. 
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(l4) 


G _ i_ 21I = £i c i c 

2A c av 2A C l c av 

The loading is known to he zero at t) = 0 ^ 1 is 

determined at three intermediate span stations. Values of 
loading at other span stations can he obtained from a load- 
ing function derived in Appendix B or, with equations (B23) 
or (B24) of Appendix B , the loading can he found at span 
positions q = 0.981, 0.831, 0.55 6, and 0.195. 


3. Rolling moment : The damping— in-roll derivative for the solu- 

tions of equations (12) or (13) is derived in Appendix B and 
given by 


PC*. 


c av 16 


0 S)t 


(is + 0.707 (Gi + G 3 ) 


(35) 


4. Spanwise center of pressure ; The equation giving center of 
pressure on the wing semi span is shown in Appendix B to he 


PCj/k 


’Ic.p. pA 


av 


~ (0.163G X + 0.248Ga + O.43OG3) 


5. 


0.082 ° lC ^ + 0 . 124 ( ^ + 0.215 f 

\ c l c av/i V G l c avyfe \ c l c av/3 

Induced drag ; The induced drag is derived in Appendix B 
and given hy 

^ = |(l^r)[ + a * s + - t g =(°i + <*>] ( 


Aerodynamic characteristics due to aileron deflection .— The solur- 
tions for the aerodynamic effects due to ailerons will hp found for three 
different spans of outboard and inboard ailerons. Cross plots of these 
data provide curves for arbitrary aileron spans. 

1. Simultaneous loading equations : The p VIl values are 

obtained from figure 1 and table I with values of H v 
given by equations (3) or (5). 
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(a) Deflected outboard ailerons : The aileron spans measured 

from the wing tip inboard are given by q a . The simultaneous 
solution for antisymmetric spanwise loading due to deflection 
of any of the three following aileron spans can be obtained 
from the appropriate set of the following equations: 


- P 11 G 1 + P 12 G 2 .+ P 13 S 3 
= P21§1 + P22<J2 + P 2 3 5 3 

= P 31 Gi + P 32 G2 + p 33 g 3 

(18) 

where G n = G n /8. 

(b) Deflected inboard ailerons : The aileron spans measured 

from the wing midspan outboard are given by i} a . The 
simultaneous solution for antisymmetric spanwise loading 
due to deflection of any of the three following aileron 
spans can be obtained .from the appropriate set of the follow- 
ing equations: 


~ PilGl + P12G2 + Pi3 G 3 
= P21G1 + P22§2 + P 23^3 
= P31G1 + P 32 G2 + P 33 G 3 


Case 

IV 

V 

VI 

n a 

0.556 

0.831 

1.000 

a l 

X 

o.o44 

0.013 

1.016 

a 2 

T 

-.017 

.961 

.9 79 

03 

T 

1.087 

1.095 

1.101 


where G_ = G n /&. 


Case 

I 

11 

III 

^a 

0.169 

0.444 

0.805 

Oi 

8 

1.003 

0.971 

0.998 

0^2 

X 

t— 

H 

O 

• 

.996 

.991 

a® 

8 

.006 

.014 

.978 


( 19 )' 
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2. Loading distribution : The spanwise loading distributions 

due to various aileron configurations include: 

(a) Full-wing— chord ailerons : The loading is known to be 

zero at tj = 0 and 1 , and is determined at three inter- 
mediate span stations. With equation (C13) and tables C 6 , 

Bl, and C7, the loading can be found at span stations 

T] = O. 98 I, O. 83 I, 0.556, and 0.195 for each of the aileron 
spans considered. With these given points and the knowledge 
that the slope of the loading— distribution curve is theo- 
retically infinite at the point of angle-of-attack discon- 
tinuity (aileron spanwise end), the loading distribution can 
be faired. 

(b) Constant fraction of wing— chord ailerons : The spanwise 

loading of constant fraction of wing-chord ailerons is equal 
to the product of the loading due to full-wing-chord ailerons 
and the effective change of angle of attack with aileron 
angle , 7 da/d5. The factor da/d5 is a function of the 
ratio of aileron chord to wing chord t = c a /c. The change 
of section angle of attack with aileron angle da/d 6 is 
presented in figure 3 , which is reproduced from figure 18 of 
reference 7 * 

Although figure 3 taken from reference 7 limits the Mach 
number range to Mach numbers less than 0.2, this limitation 
is believed to be unwarranted since theory indicates that 
da/d5 is unaffected by compressibility for the two- 
dimensional wing. However, as indicated in reference 6 , 
da/d 5 is strongly affected by low aspect ratio and will 
change appreciably if the parameter 3A becomes much less 
than twoj hence, the values of da/d5 from figure 3 appear 
to be valid for pA >2. 

(c) Arbitrary spanwise distribution o'f aileron chord : The 

aileron can be divided into several spans with constant 
da/d&, then the total loading is the sum of the products of 
the full— wing-chord loading of each span and its respective 
da/d5. 


•In using da/dS here, it should be noted that the assumption is made 
that the effective airfoil section is taken as being parallel to the 
plane of symmetry and that the section approaches a two-dimensional 
section. The validity of this assumption can be questioned! however, 
limited checks with experiment show it to be at least approximately 
correct. 
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3. B pin rip moment : The rolling moment can he found for the 

following aileron configurations: 

(a) Full -wing— chord-ailerons : The spanwise loading due to 

aileron deflection cannot he integrated with sufficient 
accuracy with equation (15)- In Appendix C, a si mil a r 
integration formula is developed that applies to each given 
aileron span. Equation (CIO) and table C5 give 


^6 

*av 


(£) 


(' 


h^Gi + + bgPa 


) 


(20) 


where for each of the cases of equations (18) and (19) the 
1^ values sure given hy 


Case 

I 

II 

III 

IV 

V 

VI 

hi 

0.l40 

0.139 

0.138 

0.146 

0.l4l 

0.140 

h 2 

.199 

.196 

.196 

.200 

.197 

.198 

h a 

.1*1-5 

.139 

.138 

.140 

.139 

.l4o 


(h) fiaastant fraction of wing-chord ailerons : For constant 

fraction of wing-chord ailerons with aileron angle measured 
parallel to the plane of symmetry, the aileron effectiveness 
is given hy 


. to 

K av d& \ K a v J 


( 21 ) 


(c) Arbitrary spanwise distribution of aileron chord : The 

deflection of ailerons for which t varies spanwise on the 
wing can he considered as an equivalent wing— twist distribu- 
tion. The effective antisymmetric twist of the wing is 
given hy 


_ _ da R 

0.41 — ' Ox 

v d6 t 


( 22 ) 
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where da/clB is now a function of spanwise position. The 
antisymmetric angle— of-attack distribution given by equation 
(22) can be divided into spanwise steps of constant angle of 
attack and the total rolling moment, can be found by the 
summation of the rolling moment due to each spanwise step. 

The rolling moments of the spanwise steps are obtained from 
a curve of rolling-^noment coefficient as a 

function of unit antisymmetric angle of attack from the wing 
root outboard. This step method is the procedure used in 
reference 1. 

A curve of / K av as a function of unit antisymmetric 

angle of attack from the wing root outboard can be obtained 
from the solutions of equation (19) for the cases TV , Y, and 
VI. An additional point can be obtained from the solution 
of case III of equation (l8), applying the relations (dis- 
cussed later) existing between inboard and outboard ailerons. 

The rolling moment due to the twist given by equation (22) 
can be obtained, by a method other than the step method, from 
the integral given by 


K av Jo d & dT l 


(23) 


which can be integrated numerically by taking the graphical 
slopes of |3C^/K av which is a function of extent of unit 

antisymmetric angle of attack from the wing root outboard. 


4. Spanwise center of pressure and induced drag ; Spanwise center 
of pressure and induced-drag integration formulas for loading 
due to ailerons are not given; however, equations ( 16 ) and 
(17) can give approximate integrations of the loading to 
obtain center of pressure and induced drag. 

5. Additional considerations : 

(a) Relation between aerodynamic characteristics for outboard 
and inboard ailerons: The spanwise loading distributions due 

to outboard and inboard ailerons bear a simple relation to 
each other. Since loading is linearly proportional to angle 
of attack, loadings are directly additive. Then, for 
outboard and inboard ailerons with the spanwise ends of the 
ailerons at the same span station. 
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G inboard • “ Tj a = l) ^outboard 


= Ci 


- Ci 


^inboard ^( t 1p = 1 ) ^outboard 


> (24) 


‘a . 


Inboard 


= 1 - 1 ! 


'a 


outboard 


These relations do not apply for Tl c .p. and Cp^ since 
these characteristics are not linearly proportional to 
loading . 

(b) Differential aileron angles : The effect of a differen- 

tial between aileron angles can be taken into account by 
considering the C of each wing panel as one— half the 
antisymmetric results of equations (20), (21), or (23). The 
total wing rolling moment is then the sum of the products of 
C ^ /2 given by equations (20), (21), or (23) and the angle 
of deflection of each aileron. Although the total rolling 
moment can be found by this procedure, the spanwise loading 
distribution can be found only approximately by the products 
of the antisymmetric unit loading G/b and the deflection 
of each aileron. However, the loading distribution so found 
will be quite accurate since this procedure neglects only 
the small change due to the induced effects of the differen- 
tially different opposite wing panels. 

(c) Aileron angles measured perpendicular to the hinge line: 
The relationship between aileron angle measured perpendicular 
to the aileron hinge line and that measured parallel to the 
plane of symmetry is given by 


tan 5 = 


tan S 
cos 


(25) 


where 

A^ sweep angle of the aileron hinge line 

5 angle measured perpendicular to the hinge line 

For constant fraction of wing-chord ailerons on straight- 
tapered wings y is given by 
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tan At = tan Ac / 4 - (26) 

where t is the fraction of wing-chord aileron measured 
from the wing trailing edge. 

Aerodynamic characteristics due to sideslip of wings with dihedral .— 
The total antisymmetric loading due to sideslip can he considered as the 
sum of that due to dihedral angle and that due to zero dihedral engle . 

For the unswept wing, the rolling moment due to sideslip for zero 
dihedral angle is generally considered negligible; however, for the 
swept wing, this effect can be appreciable. In the present report, only 
that part due to dihedral angle will be considered for the swept and 
nonswept wings . 

1. Simultaneous loading equations : The p vn values are 

obtained from figure 1 and table I with values of H v given 
by equations (3) or (5)- 

The simultaneous equations resulting from the substitution 
of 5 = Pf (see equation (ll)) and G = G/PT in equations 
(18) and (19) are applicable in the determination of the 
effects of unit outboard or inboard dihedral angle over the 
span of the ailerons considered. 

2. Bolling moment : The rolling moment due to various dihedral 

angle distributions include: 

(a) Constant spanwise dihedral angle: For dihedral angle 

constant for the entire wing semispan, the loading is given 
by the solution of case VI in equation (19) for G = G/pr 
and the rolling moment from equation (20) becomes 

— -4 = (O.l^OGi + O.I98G2 + O.lkOGj (27) 

K av* K av J 


(b) Gulled _wing : For the gulled wing, solutions of equation 

(19) for G = G/pr gives the loading, and the rolling 
moment from equation (20) becomes 

pc z - 

k — I, = (hiGi + h 2 G 2 + h 3 G 3 ) (28) 

av 1 av 

A plot of the results of cases TV, V, and VI gives the 
extent of unit dihedral angle from the wing root outboard. 
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Then, for a gulled wing, the total rolling moment equals the 
sum of products of dihedral angle of each span section «nd 
the rolling -moment contribution of the respective span 
sections. 

(c) Variable spanwise dihedral angle: If r varies spanwise, 

the rolling moment can be obtained by integration as in 
equation (23). The integral becomes 


T 1 , , a<l>VK av r) 

/ r (i> ' \ *1 

K av J 0 dr) 


d(sc lg ./K aT r) . 

where ^ r- is the slope of the curve described in 

part (b) above. 


Solution for Straight-Tapered Wings 


Charts of aerodynamic characteristics for straight— tapered wings can 
be presented in terms of geometric, compressibility, and average section 
lift-curve— slope parameters. These charts provide a ready means of 
obtaining data directly. 

Aerodynami c characteristics due to rolling .— The application of 
equation (12) for a constant value of section lift— curve slope 8 provides 
the spanwise loadings at span stations 0.383* 0.707* and 0.924 which are 
presented in figure 4 for a wide range of plan forms . . The interpolation 
formula of equation (B24) will give values of loading due to rolling at 
span stations other than those presented. With equation (15),the 
damping— in— roll coefficients 0C2pA av can be obtained and are presented 
in figure 5 for a wide range of plan forms. 

throughout the figures, K av is the constant spanwise section lift- 
curve slope or the average of a small variation. For large spanwise 
variations of k that follow the function given by equation (Bll) 
developed in Appendix B, the parameters PA/k qv and X can be replaced 

by the parameters 24 

( X 

spanwise variations of k that do not follow the curve of equation 
(Bll), the simultaneous equations for the general solution can be 
solved for arbitrary distributions of k. The H v values can be 
obtained from figure 2. 


)/(l+X) 


K m 

and —1 X, respectively. For large 
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Aerodynamic characteristics due to aileron deflection .— The appli— 
cation of equation (19), case III of equation (18), and equation (20) 
provide aileron effectiveness in the coefficient form 3Cjg/K av for 
several aileron spans. In figure 6, pC 2 gA av is plotted against extent 
of unit antisymmetric angle of attack from the wing semispan root, 
outboard for a range of wing parameters. 

As presented, figure 6 gives directly the effectiveness of full— 
wing-chord inboard ailerons for aileron spans measured from the plane of 
symmetry outboard. The effectiveness of full-wing-chord outboard ailer- 
ons for aileron spans measured from the wing tip inboard is given by 
figure 6 directly by the relations of equation (24). For full-ving- 
chord ailerons located arbitrarily on the wing semispan, the aileron 
effectiveness' can be obtained directly from figure 6 as indicated in the 
following example sketch. 



With the full-wing— chord values given above, the effectiveness of 
constant fraction of wing-chord ailerons or ailerons of arbitrary spanwise 
chord distribution can be found through use of equations (21) or ( 23 ) with 
the da/d5 values of figure 3. 

Aerodynamic characteristics due to sideslip o f wings with dihedral. - 
The application of equation (19), case III of equation (l8), but with 
5 = pr, and G = G/pT, and the use of equation (28) provides rolling 
moments due to dihedral angle for the wing in sideslip. These rolling 
moments are given in the coefficient form PC 2 j 3 A av r which is the s ame 
function of q as $C l /K av and is presented with pC2 6 A av in figure 6. 
Figure 6 with equation (29) will provide the rolling moment due to 
sideslip for any symmetric /spanwise distribution of dihedral angle. 

For dihedral angle constant spanwise , the rolling moment is given by 
the value at q = 1 in figure 6. These values for constant spanwise 
dihedral angle are presented in figure 7 as a function of aspect ratio 
for various values of sweep angle and taper ratio. 
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DISCUSSION 


Effects of plan— form parameters on aerodynamic characteristics for 
straight-tapered wings are shown by plots against the various parameters. 
Compressibility is discussed and formulas given for a range of plan forms 
at sonic speeds. Theoretical considerations and experimental comparisons 
indicate the order of reliability of the present theoretical results. 


Straight-Tapered Wings 


The spanwise loading distribution due to rolling for several plan 
forms is presented in figure 8. These curves are the result of applying 
figure 4 and the loading interpolation formula of Appendix B. The 


loading coefficient is given as 


^c.p. f c l c \ 
(''lc.p.^-Q V c Z c av^ 


to make the total 


loading on the semispan constant and thus show more clearly the changes 
of distribution due to sweep and taper ratio. Figure 8 shows large 
c hang es in loading distribution for the zero tapered wing. The effects 
of sweep are generally as expected, namely, that sweepback shifts the 
loading outboard . 

Effects of plan form on the rolling moment due to rolling is shown 
from cross plots of figure 5 which are presented in figures 9 and 10. 
For higher aspect ratio, figures 4, 9, and 10 show the marked lowering 
of rolling moment due to sweep. Figure 9 indicates that for low aspect 
ratio, the rolling moment becomes essentially independent of sweep and 
taper. The taper effects on rolling moment as seen in figure 10 are 
small except for values of taper ratio less than 0.25. 


Typical spanwise loading distributions due to full-wing— chord 
aileron deflection are shown in figure 11. These curves were faired 
with the aid of the loading interpolation function of Appendix C and, 
at the aileron spanwise end, care was taken to make the slope large. 


Wing geometry effects on aileron effectiveness for full— chord 
outboard partial— span ailerons (with aileron angle measured parallel to 
the plane of symmetry) are given in figure 12. The geometry effects on 
PCjg/Kav are similar to those on the damping— in— roll coefficient. Com- 
parison of figure 12(a) with figure 9 shows that Cjg approaches the 
zero— aspect— ratio value in the' same manner as does C^. Figure 13 gives 
comparative effectiveness of inboard and outboard ailerons for swept 
wings. As sweep increases, the difference of effectiveness between 
inboard and outboard ailerons decreases showing that inboard ailerons 
for highly swept— back wings approach the effectiveness of outboard 
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ailerons. Since da/dS 1560011168 large rapidly at small values of t 
(fig. 3), then, for a given aileron area, narrow full— span ailerons for 
swept-back wings may be more desirable than larger-chord outboard 
ailerons. The relative effects of figures 12 and 13 apply equally well 
for constant fraction of chord ailerons, since the data would differ only 
by a constant factor da/d&. 


Compressibility 


From the three-dimensional linearized— compressible— flow equation, 
it can be shown that the effects of compressibility will be properly 
taken into account if the longitudinal components of a wing plan form are 
increased by the factor 1/3 . Or, alternatively, if the linearized com- 
pressible flow equation be divided through by P 2 , then the lateral and 
vertical components of a plan form are decreased by the factor p. In 
both cases, the incompressible local lift is increased by the factor 1/p 
and the compressible local lift coefficient can be written as the 
parameter PCj. 


With these relations known, an incompressible theory can be made 
into a compressible theory subject to the limitations of the linearized 
compressible flow equation. The geometric parameters of a wing are 


simply substituted by PA,Ap = tan -1 - A and Pb. With local lift 
coefficient given by Pcj, the dimensionless loading becomes 


3°2C = cjc 
2pb 2b' 


The wing-chord distribution remains unaltered. 


The sonic speed results of reference 6 can be used as a limit point 
in the present theory for a curve of the variation of antisymmetric 
aerodynamic characteristics with Mach number. The following equations 
apply at the speed of sound to plan forms with all points of the 
trailing edge at or behind the upstream line of maximum wing span: 


For outboard ailerons. 



— flA 
32 


6 

For inboard ailerons, 

C 


Cift = £ sin 3 0, where = 1 — cos 9 


'lb 


= 4 ( 1 — sin 3 0), where T) = cos 6 
o a 
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Eeference 6 shows that aileron effectiveness at the speed of sound is 
independent, of the chordwise location of the aileron hinge line, pro- 
vided the hinge line remains ahead of all points of the trailing edge. 


Accuracy of the Seven— Point Solution for Ailerons 


The prediction of aileron effectiveness for given aileron spans 
with wing twist determined by zero-aspect— ratio theory at only seven span 
points to satisfy the boundary conditions has been theoretically shown to 
be sufficient by comparing results with the computation of a typical 3.5 
aspect ratio, 45° swept wing with 15 span points satisfying the boundary 
conditions. The process of finding aileron spans for the 15-point method 
was the same as that in Appendix C. The curves showing the variation of 
C z 6 ..with aileron span for the 7- and 15-point computations were identical. 

• 

The solution for the angle-of-attack distribution that includes a dis- 
continuity can be compared with the solution for the continuous angle-of- 
attack distribution by considering an aileron such that the angle-of- 
attack distribution is equivalent to that of the rolling wing. The damp— 
ing-in-roll coefficient then can be found by use of equation (23) which 
reduces to the form 


n 1 dCjg 

for a = — J , and integrating by parts 

Cip =jf Cl a d V- C Vi 

This relation states that Cjp is equal to the area between a curve of 
figure 6 and the line of Cjg for tj = 1. The curves of figure 6 were 
found by the simplified lifting— surface theory with antisymmetric twist 
determined by zero-aspect-ratio theory. The values of Cj_ obtained in 
this ma nn er from figure 6 were identical to the Cj_ values given by 
simplified lifting— surface theory for continuous linear antisymmetric- 
twist distribution. 

As further theoretical check, the values of rolling moment due to 
constant spanwise dihedral angle are obtained from 15— point computations 
in reference 5 for taper ratio equal to one, with which the present 
theory for the 7— point method is in exact agreement. 
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Comparison of Theoretical and Experimental Results 


The electro-magnetic analogy method of reference 8 provides damping- 
in— roll coefficients for an aspect— ratio range of unswept , tapered wings. 
The results of the present theory and those of reference 8 are compared 
in figure 14. Except for the taper ratio effects on Ci p the comparison 
is good. The rounded-wing-tip values of Cj p given by NACA Rep. 635 
(reference l) are included in figure 14. Since rounded wing tips gener- 
ally give values of Cj p about 6 percent lower than straight wing tips, 
the values of MCA Rep. 635 appear to be appreciably too high for lower— 
..aspect— ratio wings. The present theory and the theory of reference 8 
approach the value given by the zero-aspect-ratio theory of reference 6 
quite satisfactorily. The results of the present theory may be further 
assessed by the comparison with the results of low-speed experiment as 
given in figure 15 for the range of plan forms presented. For further 
experimental verification of the accuracy with which Cj can be deter- 
mined by the present theory, the reader is referred to reference 5 which 
supports the theory as well or better than figure 15 of the present report 

The loading distributions due to rolling as given by the present 
theory are compared in figure 16 with low— speed experimental results for 
a range of swept wings. The sweep angle seems to have considerably more 
influence on loading distribution as given by experiment than the theory 
indicates. The experimental pressure data, however, were very erratic 
and no firm conclusion can be made. 

Experimental values of rolling effectiveness due to aileron deflec- 
tion are compared with theoretically predicted values in a correlation 
diagram given by figure 17. Included are the results of a wide range of 
plan forms which do not vary consistently with any geometric parameter or 
aileron configuration. Sketches of the plan forms and ailerons are drawn 
about the points of correlation. The theory makes use of the curve of 
figure 3 giving da/dS for a sealed— gap aileron over a range of deflection 
of ±10°. Experimental results for aileron deflections greater than 15° 
measured perpendicular to the hinge line were not included. The correla- 
tion points of figure 17 scatter appreciably; however, the mean line of 
the points does approximate the line of perfect correlation. 

Figure 17 does not account for effective plan— form change due to 
da/d6. Only the effectiveness of the low— aspect— ratio triangular wing of 
figure 17 is exceedingly in error, which is the result of neglecting plan- 
form change. 

The plan-form change due to da/d5 can, in part, be considered anal- 
ogous to that due to section lift-curve— slope change. Thus, the total 
section lift of a wing,, the chord of which is reduced by da/d5 and which 
is at an angle of attack 6, is equal to the. lift of the wing-aileron 
section for which the aileron only is deflected at the angle 6. This 
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change in plan form, unlike the section lift-curve— slope change for which 
the chordwise loading remains constant, does not account for a large 
change in chordwise loading. If the lifting line is considered to he at 
the chordwise center of pressure, then, for partial-wing-span ailerons, 
the lifting line is in effect broken at the aileron spanwise end and the 
present theory becomes invalid. For the case of full-wing-span ailerons, 
the lifting line in effect remains unbroken and lies along the center of 
chordwise pressure. For this case the wing chord can be reduced by da/d6 
to account for plan— form change; however, although in the limit of zero 
aspect ratio the results are the same as those of reference 6, this pro- 
cedure does not with sufficient accuracy account for the chordwise loading 
shifting aft at intermediate aspect ratios. For control surfaces, the 
effective plan-form change due to da/d6 is appreciable for the low- 
aspect— ratio wings such that in the limit of zero aspect ratio the span- 
wise loading is independent of the ratio of aileron chord to wing chord 
(reference 6). However, for moderate aspect ratios, da/d6 can be used 
without accounting for plan— form changes as comparison with experiment 
indicates. 

Experimental values of Cjjj/r are not compared with the present * 
theory since reference 5 gives ample support of the theory. 

CONCLUDING REMARKS 


The determination of antisymmetric loading for arbitrary wings is 
shown to be easily obtained by the solution of three simultaneous equa- 
tions. The coefficients of the simultaneous equations are presented in 
charts of parameters that include wing geometry, compressibility, and 
section lift-curve slope as arbitrary quantities. Thus the loading for 
an arbitrary antisymmetric angle-of-attack distribution can be Bimply 
found once the angle— of— attack distribution is chosen. 

For the important cases of antisymmetric loading, roll, and aileron 
deflection, the angle-of-ettack distribution is given and the simultaneous 
equations are formed. Loading for these cases can be found by simply 
obtaining from charts the coefficients corresponding to the wing geometry, 
Mach number, and lift-curve slope, inserting in the appropriate equations 
and solving. 

Integration formulas for the loading distributions are given which 
enable the aerodynamic coefficients and Cjg to be found. The 

rolling moment due to sideslip of a wing with dihedral is shown to be 
equivalent to that of aileron deflection and a procedure for determining 
its value is given. 

For the special, case of straight— tapered wings, the loading distri- 
butions and values of and Cjg are given in the chart form for a 

range of wing plan forms. 
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Experimented, and theoretical verification of the theory is shown to 
be good. The theory is applicable for large aerodynamic angles, provided 
the flow remains unseparated. The compressibility considerations are 
reliable to the speed of sound subject to the limitations of the linear- 
ized compressible— flow equation. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., Dec. 22, 19^9. 


APPENDIX A 

EQUATIONS FOE THE DETERMINATION OF ANTISYMMETRIC LOADING 


Unsymmetric Loading 


From NACA TN 1V76 (reference 2), the aerodynamic loading is obtained 
by solution of the linear simultaneous equations 


% *( 0 „ - 


n G v 


m 1 

^**vnGn, ^ = I* 2, . . . m 

n=l 


(The prime indicates the value for n=v is not summed.) 


where 

It 


G = 


bV 


b* v = 2b vv + 


b* = 2b 
Vn Vo. 


(A)' 

0 ) 


VV 


g vn 


(Al) 


(A2) 

(A3) 


b vv and b-^ are coefficients independent of plan form. 
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g, 


— 1 r L ( V '°) f no + L(V,M+1) f. 


Vn 2(M+1) 




n,M+i 


d=l 


u)f. 


up 
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2 


coefficients independent of plan form. 
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The above equations involve computations over the entire' wing. 
However, if the loading is assumed to be symmetric or antisymmetric, the 
computations can.be reduced to less than half the work. The case of 
symmetric loading is developed in reference 2 and the antisymmetric case 
is developed in the following section. 


Antisymmetric Loading 


For antisymmetric loading, the loading on each side of the wing has 
the same magnitude and distribution but with opposite sign, or 


a V “ ““m+i-v 


G V G m+i-v* or G n G m+i-n 


(A6) 


Equation (Al) can then be written as 


m-1 

2 


“v - < b *v -»Vi- v )°v b *vn- b *v, a+i-n) 

n=l 


(A7) 


where the summation is only to since G m+1 = 0 for antisymmetric 

loading. 2 

With equations (A2) and (A3), equation (A7) becomes 


Oy = 


2 (b vv -b v , m+1 _v) +(|^)(Svv _ 8v,m+i-v) 


G\> “ 


m-1 


* 


n=l 


^ ^ ” ^cy^ (^®Vn ®V,m+i-n ^ ^n 


(A8) 
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Now, from equation (a4) 


8vn“Sv,m+i-n “ 


-1 


2(M+1) L 


fl( V ,0) (fno-W-n.o) + 


L( V,M+l) (f n<M+1 

M+i^ 


where 


M 

Y ( f : 


np""^m+i-n, p) 


P=1 
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DP 


= Pi si.n pa<Pn cos pi<Pp 

Pi=l 


_ nit an 

Cp = y CPiI = 

n m+l ^ M+l 


From equation (A9) , f^— : f m+i-n,p can defined as 


f * = f -f 

np np m+i— n , p 


then, using equation (A10) 


f* 


W ~ 55T 


m 

Y Pi cos Pi^ (sin mcp n - sin Pi^m+i-n) 


Pi=l 


= Y ^ cos ^ lCp P sin ^ lCp n ( 1+ cos Pi*) 

Pi=l 


and, since the terms of the summation for odd pi vanish, 
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f V m+l 


I 
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Pi sin nxtp n cos picp^ 
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From equation (All) , 

f*nn = **n,M+i-n (A12) 


Combining equation (A9) with (A12) and defining 


®*Vn ^Vn~®v,m+i-n 


then 


g*vn = ~r— 
2(M+1) 


M+l 

2 


I 

|j.=0 


L(v,p) + L( v,M+l-p) 



(A13) 


where for jj. = 0 and >f*^ e Q. ua f to half the values given by 

equation (All) in order that the products can he fitted into the summa- 
tion. With equation (A13), equation (a 8) can he written as 


m-1 

2 , 

®V = (2C V + ^ g*vv) G v — ^ ( 2C Vn “ S*vn)°n (Al4) 

n=l 

— ^v^m+i-v 

— ^Vjm+i-n 

From reference 2, 

b sin cp n f l-(-l) n ~ V ~ 

Vn (cos 9 n - cos cp v ) 2 L 2(m+l) _ 

which gives zero values for b Vn for even (n— v) values. Then, since 
m is odd, 


where 


V = 1,2,3,... ^ 


C V — 


C Vn - b vn 


h =0 

V, m+i— v 


C v = b 


W 


and 


32 


NACA TN 2140 


It should he noted that L( v,p) simplifies somewhat, for the anti- 
symmetrically loaded wing since T) now is only positive in equation 
(A5). If only positive values of Tf sire used, then equation (A5) can 
be written as 

L*( V,n) = L* = L(q,q) + L(t),-^) 

= L(v, m ) + L(v,M+l— p) 

In summary, the foregoing analysis for the antisymmetrically loaded 

wing gives , 

ffirrl 

2 

I Pvn^n (A15) 

n=l 

V= 1,2,3,... 


2b vv + T" g *vv for n = v (Al6) 
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V C J X T, v“p ' 


)(TUrV tan A] 2 +(k-) 2 ( Tlv ^) 2 


1 f / [1+( I^.K^P )tan ~ (^) 2 ( yV‘ 




1 + 2 (E_) T| tan A 
0 v v 


-4 


- 1 


2 tan A 


J [i+(cv: 


cv^v tan A] + ( C y) 


l+2(— )Tj v tan A 

C -VI 


T) = cos 9 V where cp v = 
v ^ m+1 


ti = cos cp where cp 
V H M- 


M+1 


Tiv, = cos cd where cp = 

' n • % n m+1 


For a discussion of the relative accuracies obtained for a choice 
of values -of M and m, see reference 9. The most favorable applica- 
tion is with M=m. 
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APPENDIX B 

DERIVATION OF RELATIONS USED IN THE METHOD 


Application of Appendix A 


With Appendix A, the anti symmetrical loading on a plan form for any 
anti symmetrical distribution of ay can be found. The principal work 
in the computations is to obtain the coefficients of the simultaneous 
equations (A15). These coefficients can be presented in charts for the 
complete range of geometric plan— form parameters into which are intro- 
duced the effects of compressibility and section lift-curve slope. With 
the loading due to rolling known, the coefficients and derivatives are 
obtained by integration formulas. 


Section lift-curve-slone effect .- For a two-dimensional wing with 
the loaded line at the quarter— chord position, the position x aft of 
the loaded line where the induced downwash equals the angle of attack of 
the wing can be obtained by the Biot Savart Law as 


or 

or 


T r c j cV 

w = where T c = — — 

2nx 2 

W C Z C q 

v 4-itx 
dc z _ W 

da c 


then 


* - <H P- 

da 


where dcj/da is the section lift-curve slope. Two-dimensional section 
compressibility effects that do not follow the Prandtl-Glauert rule can be 
given consideration by taking the ratio of Uc z /da) compressible at a given 

Mach number to 2it/p. Let k be the ratio of the section lift — curve slope 
at a given Mach number to 2n/p or (dc z /da) compressible = 2«n/p, then 


x = k(c/2) ' 

Then the induced angle, k (c/2) aft of the loaded line, is equal to the 
angle of attack of the wing. For K = 1, this is at the three-quarter- 
chord line. For section lift-curve slope less than 2n, k is less t han 
one and the downwash is equal to the angle of attack at some point between 
the one-quarter— and three-quarter-chord line. 
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To take into account the section lift-curve— slope variation in the 
present theory, the downwash must he found at a distance k(c/2) aft of 
the xojg.de d line. From the formulas of the summation in Appendix A, 
b/c v 'should be taken as b/KyCy, where k v is the ratio of section 
lift-curve— slope for a given Mach number at span station v, to 2n/3. 

Derivation of parameters for Pvq* - The Pvn coefficients, as 

defined by equation (Al6) in Appendix A, depend on plan— form geometry in 
the (b/cv)!*^ functions only, or p^ is a function of b/c v and 
sweep angle. As previously shown, b/c v is also a function of the span- 
wise variation of section lift-curve slope and is effectively equivalent 
to b/K v c v , where k v is the ratio of section lift— curve slope for a 
given Mach number at span station v to 2 ji /| 3. The p^ coefficients 
can be plotted against b/K v Cy with sweep angle as a parameter; however, 
b/ K v c v will vary from zero to very large values for a range of plan- 
form geometry, and the plots become unwieldy. For a range of aspect 
ratio, the values of b/ityCy are a maximum for the zero tapered wings 
when tj v >0.5 (provided plan— form edges are not concave) and a maxi mum 
for the inversed-tapered wings for tj v < 0.5. The ratio of b/KyC v 
for q v > 0.5 for any plan form to those of the zero tapered wing or 
the ratio of b/Kyoy for v < 0.5 for any plan form to those of the 
inversed-tapered wing gives a geometric parameter for any plan form that 
has maximum values that depend only on aspect ratio. 

The chord distribution for straight— tapered wings is given by 


b a(1+A) 

cy = 2 [1- fo v |(l— A) ] 

Then, for A = 0, 


(Bl) 


b _ 1 

ACy 2(1- h v |) 

and for A = 1.5 

-B = 5 

Ac v 2(2+ |n v | ) 


(B2) 


(B3) 


The ratio of b/nyc v 
a geometric parameter as 


to equations (B2) and (B3) gives, respectively. 


b/lCyCy 

(b/Acy)^ 

b/kyCy 

(b/Ac v )^ =1> 5 


= 2(l+ lv )^ 1 ^)for 0.5 <T|y < 1 


2( 2+r| y ) / h 


(^r v ) for °s"vso. 5 _ 


(b4) 


5 
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Let Hy be defined as two-fifths times- the values of equation (B4) 
(the fraction two— fifths is introduced to give Hy the approximate 
values of Py_ to simplify plotting procedures ), then adding effects of 
compressibility (see Discussion section) 


Hy 



(B5) 


where 


d v = 


Mi-n v ) 

5. 


for 0.5< 


<1 


M2+t] v ) 

25 


for 0 <rj v < 


0.5 


For tapered wings, Hv simplifies to 


H v 


2(l-^i v )(l+X) 

5[i-n v (i-x) ] 



for 0.5 —T v < 1 


2(2+t) v )(H-\) 

25 [1 — T1 V ( 1—1 ) ] 



for 0 <tj v <0.5 


(B 6) 


j 


Plots of pyjj against H v in the range of H v = 0 to 4 will give 
Pyn coefficients for wings of any chord distribution for aspect ratios 
up to 10 or 12. 

Linear asymptotes of p vn .- For large values of H y , the p yn 
functions become linearly proportional to Hy. Since this linear char- 
acteristic appears at relatively low values of H v , the simply linear 
relation between p Vn and H y is quite usable. 
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The L* Vn function of Appendix A is multiplied by b/cv and the 
product is linearized. 






With the values of equation (B7) substituted into equation (AL6) 
the values of for arbitrary sweep angle are obtained. Thus, 

for m=7, the following equation (b 8) gives values for p vn as 
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•ll 


■12 


( 3, i — ■ + 1.026 tan AJ Hi + 7.968 - 1.494 sin A + 0.014 tan £ + 
\cos A / 2 

0.08e f ^■Z 1 . + 0 -i»l<? , jgsA) - 0.068 fe ^ l±°- W7 tag 2 A \ 

. \ tan A / \ tan A / 

0.034 ( W . 1 . ^ 

\ tan A / 

( - • ■ 8 — - 2.901 tan a) Hi - 3.138 + 1.080 sin A - 0.034 tan § - 
\ cos A / ^ 

1-a/i + 0.0016 tan 2 A ^ _ 0 096 ^ 


tan A 


1- 4/1+ 0.1717 tan 2 A 
tan A 

A 


•13 


( -0.1] 6_ + 1>026 tan A j Hi + 0>129 _ 0>869 sinA + 0<o82 tan A + 
\cos A ' 2 

0.014 A-V 1 + 0.0016 WAV 0 068 A- + 0.0177 tonA. 

\ tan A / ' tan A 

0.034 A / 1 + °- 1T17 jg l) 

V tan A ' 


/ 0.221 

21 \cos A 


+ 0.534 tan A 


0.088 f - 


J 1 


a)h 2 - 2.088 - 0.383 sin A- 0.018 tan | - 
) - 0.044 (tJ± + °‘ 0294 tan2A ' 


+ 0.0177 tan A 


tan A 


tan A 


0.037 


1— yl + 0.0886 tan 2 A 
tan A 


) 


^22 


2i212 h 2 + 4.596 - 0.146 sin At-0.o 44 tan £ + 
cos A ^ 

0.125 (— - 1 + °- 017T tan2A N ) o. 0 i|i| Azj/TT 0.0294 

\ tan A / V tan A 


tan 2 A A _ 


/l- V 1 + 0.0886 tan 2 A^ 
^ \ tan A y 
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A ). 221 „ , 

r 23 = l 7 - 0.534 tan A 

23 Vcos A 


l) H 2 - 1. 


912 + 0.221 sin A + 0.107 tan — - 

2 


0.044 (- -iL 1 + 0,01 77 tan 2 A \ t 0>0 i8^ jc. ^ 1 + °- Q2 94 tan 2 A \ + 
\ tan A / \ -tan A J 

. . / 1 — J 1 + 0.0886 tan 2 A \ . 

0.044 ( ^ s ) 

V tan A / 


‘) 


p = f ~0»028 _ o.l64 tan a) H 0 + 0.149 + 0.324 sin A + 0.034 tan - - 
3J - \cos A 


0.082 1 + 0.1717 tan 2 A \ Q>1 g, A- 7~1 + 0.0886 tan^ A\ + 

\ tan A J ' \ tan A ) 

\-j 1 + 0.1993 tan 2 A \ 


0.197 


tan A 


V 


= (O-ffi- + 0.464 tan A ) H a - 1.570 - O.389 sin A -0.082 tan ^ 
32 Vcos A / 3 2 


0.231 


l-J 1 + 0 . 


1717 tan 2 A 


tan A 


^-0.082 


1- yrr 0.1993 tan g A 
tan A 


) 


Poo - 0.164 tanA ) H 3 + 3.417 + O.083 sin A+ 0.197 tan 4 - 

33 Vcos A / 2 


0.082 ( - + 0 • 1T17 tan 2 A V + q ^3 f 1~ *J 1 + 0.0886 tan 2 ~ 'V + 

\ tan A J \ tan A / 

0.034 ( l + 0,1 993__tan 2 aV (b8) 

V tan A / 


Linear spanwise distribution of ( icc) v /(icc) av .— With the condition 
that the product of section lift-curve slope ana wing chord varies 
linearly spanwise, then 


KC 


2b 

At ] 


“^1“ 1~ (Krp /Kp ) X ] 


} 
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and equation (3) becomes 


h v = a v (pa k ) 


l+( Kp/ Kp)X 

2 | l-n v [1 -(k t /k r )1]| 


(B9l 


where A K is the aspect ratio based on the wing chord equal to kc. 

In equation (B9), Hv is reduced to terms of two parameters. Expressions 
of A k in terms of aspect ratio for straight-tapered wings and the dis- 
tribution of section lift-curve slope can be found. 


For straight— tapered wings 

2b 

A = 

©^(l+x,) 


nnrl since kc is linear 


2b 

Ak= VrI^vV ’ 11 

then 


cvW 1+x 


and equation (B9) becomes 


PA 

l+( Kq-i/ 

( + KijA ) / 1+\ 

2 1 

1-11 V (1-( k t/ k r)^^ 


(BIO) 


The distribution of k for straight-tapered wings is given by 
(kc) v _ 1-t) v [1 -( k t/ k r)^] 

Khj — - — ■ ■■■. r.- = “* 


K B 


i-n v (i-x) 


(Bll) 


Equation (BIO) is in terms of two parameters given by ! 

j_ (Kj^+KpX)/l+\ _ 

and ( Kp/ Kp ) X . Solutions for spanwise loading in terms of these two 
parameters and Ap are valid for the distributions of section lift- 
curve slope given by equation (Bll). Equation (Bll) indicates that at 
X = 1 , kv is a linear function and at X = 0, k v is a constant. 
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For values of X between 0 and 1, k v is a curve in the region 
between the linear function and a constant. 


Equation (BIO) is given by figure 2 for m=7> but with the ordi— 

H v 


nate given by the parameter 


( k t A e )x. 


|3a/[(k r -h< t \)/(1+X)] 


and the abscissa by 


For the case of linear distribution of (kc)v and straight- 
tapered wings for which the chord and section lift-curve slope can be 
specified in three parameters, the loading and associated aerodynamic 
characteristics can be presented for a range of the parameters Ap, 

3A/[(k r +k t X)/( 1+X)] , and (kt/k e )\. 


Integration of Antisymmetric Loading 


Rolling-moment coefficient and derivatives.— Rolling-moment coef- 


ficient is given by 


P c l = 7T J G(q)qdri 


(B12) 


where 


T| = cos 9 


which, by an integration formula, 

m 


J f(q)dq = ^ f(q n ) sin cp n 

_1 n=l 


(B13) 


becomes 


PC, = 


npA 


1 2(m+l) 


G cos <t> sin cp 


n=l 


7tpA 

4(m+l) 


m 


n=J 


G n sin 2cp n 
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Since the loading is antisymmetric, G m+1 = 0, and 


m— 1 

2 


I G “ 8111 ^ 


n=l 


(B14) 


Tor spanwise loading due to rolling, the loading is found as a 
function of pb/2V, then equation (Bl4) divided by pb/2V gives 

m— 1 
2 

I ^ sln 2<Pn ' (B15) 

n=l 

where _ 

G = G/(pb/2V) 

The rolling moment due to. ailerons will be found in Appendix C. 

Induced drag .— The induced drag coefficient is, with equation (B13), 
given by 

pi ,t|3A V 

PCp. = PA I aiGdq = VHy sin 9 V 

_1 V=1 

where a^ v is one-half the induced angle of the wing wake given by 
equation (Al4) for c v = then for antisymmetric loading 


m-1 

2 


PCDi = 


2ffpA 

m+1 


•K 

V=1 




m-1 

■t* 

n=l 


vn G n ) sin 4> v 


(B16) 


where the prime indicates the value of n = V is not s umme d - 

Spanwise center of pressure .— The center of pressure on the wing 
half panel is given by _i 


/ GqdTi 
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The numerator is equal to PC^/PA. If the Fourier series for loading is 
assumed. 


m 


c<<p) . y a sin cp, the denominator becomes 


^i=even 


V r- V — 

. Zj a m / 2sin 'P sin dlCpdCp = ) 2 

Hi=even ° M-i=eveh 


Ui 


di^-1 


then 


pc 7 


'c.p. 


f 1 — / \ 

■“'I- ^ (_ 1 ) 2 fe) 


(B17) 


Hi=even 


where , a ( j. 1 are the Fourier coefficients. 


Loading-due— to— rolling function and interpolation table .— The 
Fourier series that approximates the antisymmetric loading with only a 
few terms is given by 


G(cp) 



Hi=even 


sin (iiCp 


(B18) 


The loading G n is determined at span positions of ti = cos 'cp_ where 
nn 11 

9 n = The a^ are given by 



sin nicpdcp 


(B19! 


With the quadrature formula of equation (B13), equation (B19) becomes, 
for antisymmetric loading. 


k 

a ^ 1 - m +1 


m— 1 



n=l 


sin mcp n 


(B2( 
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For m = 1, the a (il coefficients are equal to (for even Ma) 



Equation (Bl8) with (B2l) can he arranged to give 

G (cp) = sin 2 cp + sin 4 cp + 

2 \ 2 

i (sin 2 cp — sin 6 cp )Gz + 

— ( sin 2 cp — sin 4 cp + 

2 ^ 2 


J~2 


sin 6 


<p) Gi +1 


> (B22) 


^ sin 6 cp^G a 


With equation (B22) the loading due to rolling can he determined at any 

W ' 

span position. Letting cp = cpjj. = "g and tabulating the factors of G n 

as e n k > an interpolation table may he obtained to determine loading at 
span station k. 


TABLE Bl, e ^fr 

[m=7] 


'lk 

0.981 

0.831 

0.556 

0.195 

BS1 

1/2 

3/2 

5/2 

7/2 

i 

0.8155 

0.5449 

-0.1622 

0.1084 

2 

-.2706 

.6533 

.6533 

-.2706 

3 

.1084 



.8155 


where 

n=l 


enk G n 


(B23) 
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Equation (B23) may be used for interpolation of any form of loading coef- 
ficient, thus 


( 


_c*c_\ 

^Z c av/ 

k 



(B24) 


APPENDIX C 

DETERMINATION OF ANTISYMMETRIC WING TWIST FOR FINDING 
SPANWISE LOADING DDE TO AILERON DEFLECTION 


Wing Twist for a Given Aileron Span 


The determination of loading for an angle-of-ettack distribution 
that contains a discontinuity by a method which satisfies the boundary 
conditions at a finite number of points can be made by increasing the 
number of points until the solutions become sufficiently accurate. For 
the method as given in Appendix A, the number of points that satisfy the 
boundary conditions is given by m. For the large value of m required 
for accurate results, the computations become exceedingly laborious; how- 
ever, a procedure using a moderate value of m can be determined by use 
of a low-aspect-ratio theory with which a wing twist can be found that 
duplicates the results of the discontinuous angle-of-attack distribution. 

A theoretical but relatively simple method of finding spanwise load- 
ing due to inboard and outboard ailerons for wings of low aspect ratio is 
given by reference 6. In the present theory, as aspect ratio approaches 
zero, g* Yn values of Appendix A become zero and the pvh coefficients 
given by equation (Al6) become constant or independent of plan— form shape 
and equal to 


p Vn ~ -2C Vn 

Pw = 2b vv 


(Cl)- 


These coefficients are given by the relations under equation (A15) and 
Pvn can be tabulated. 
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TABLE. Cl.— p-^ 
[For m=7 and A=0] 


n 

1 

2 

. 3 

i 

10 . 4524 

-2.0000 

0 

2 

- -3.6954 

5.6568 

- 1.5308 

3 

0 

-2.0000 

4.3296 


With equation (A15) , antisymmetric loading can be found for zero-aspect— 
ratio wings. As a comment on the accuracy of the present theory for 
m=7 , the solution of equation (A15); with the A=0 p Vn values for load- 
ing due to rolling gave the same values at the three semi span stations as 


does reference 6, namely, G(cp) 


( pb / 2 V ) sin 2cp 

4 


The zero-aspect— ratio theory of reference 6 shows that all span 
loading characteristics are independent of plan— form shape for zero 
aspect ratio. This independence makes that theory ideal for obtaining 
the boundary conditions of the present theory for zero aspect ratio, which 
should apply with the present theory for higher aspect ratios for which 
plan— form shape has an effect on spanwise loading. The boundary condi- 
tions of the present theory are. given by the antisymmetric values of <x v 
in equation (A15). The problem is to find what anti symmetrical distribu- 
tion of coy is required for the present theory to duplicate the exact 
loading distribution given by reference 6 for a given aileron span. 


The aileron spans are arbitrarily chosen for the present theory as 
the mean value of the spanwise trigonometric coordinate of the downwash 
point at a section angle of attack equal to zero. For m=7, three 
aileron spans can be defined for both outboard and inboard ailerons. 

Let q a be the aileron span, and 8 the spanwise point of the end of 
the aileron, then 


n = 1 — cos 8 for outboard ailerons 

'a 

n = cos 8 for inboard ailerons 

'a 
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For the present theory, the aileron spans defined are tabulated as fol- 
lows : 


TABLE C2 



Outboard 

Inboard 

Case 

1 

11 

hi 

TV 

V 

VI 

0 

3 £ 

5 # 

7« 

5* 

3* 



16 . 

16 

16 

16 

16 

0 

^a 

0.1685 

0.4444 

0.8049 

0.5556 

0.8315 

1.0000 


For the aileron spans listed in table C2, the exact span loading 
distribution can be found from reference 6 . With the p^ values listed 
in table Cl and the exact values of Gi, G 2 , and G 3 from reference 6 , 
equation (A15) gives the twist required for the present theory to give 
the loading distribution for each case listed in table C2 or 


r- •“;' l 52 t (r)- 3 - e 554 (r) 

r-(r>’«(rHr) 

= - 1.5308 4 . 3296 ( 52 ) 


> (C2) 


j 


The spanwise loading distribution from reference 6 for outboard ailerons 
is given by 


G(q>) 

8 


outboard 


( cos cp — cos 0) in 


( cos q> + cos 0) in 


sin 


sin 


cos 


cos 


0+9 


0-9 


0+9 


0-9 


(C3) 
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For the full-wing— span aileron. 


e = | or T} a = 1 


g(<p) 

8 "^ a = 1 


= — cos cp Zn 

it 


1 + sin cp 
cos cp 


(C4) 


For inboard ailerons, with the same value of 0 


p(cp) 1 


row] 



8 

inboard 

6 

hn=l 

6 


outboard 


(C5) 


With equations (C3), ( C4) , and (C5), the spanwise loading Gi, G 2 , 
and G 3 at span stations cp = it/8, it/4, and 3«/8, or q = 0.9239, 
0.7071> and 0.3827 can be tabulated for each of the cases given in table 
C2. 


TABLE C3 


2a 

8 

Case I 

II 

1 

III 

IV 

V 

VI 

Gi 

8 

0.1136 . 

0.1919 

0.2316 

0.0454 

0.1237 

0.2373 

g 2 

8 

.0500 

.2800 

.3851 

.1164 

.3464 

.3964 

i 

8 

.0190 

.1022 

.3620 

.2922 

.3754 

1 

.3944 


The twist distribution required for each case is obtained with 
equation (C2) and table C3, tabulating 

TABLE C4 


a 

8 

Case I 

II 

III 

IV 

V 

VI 

a-i 

8 

1.0029 

0.9713 

0.9979 

0.0444 

0.0128 

1.0157 

og 

6 

.0174 

.9957 

.9913 

-.0169 

.9614 

.9788 

03 

8 

.0056 

.0139 

.9777 

1.0868 

1.0951 

1.1007 


With the twist distribution given by table C4, equation (A15) can be used 
to solve for spanwise loading due to ailerons for any of the six cases. 
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Rolling Moment Due to Aileron Deflection 


The rolling moment is given by 

it 

G(T) sin 2 cpdtp (C6) 



For span loading due to ailerons, the loading distribution is distorted 
sufficiently such that the quadrature formula given by equation (B 13 ) is 
not sufficiently accurate for m = 7 to integrate equation (C6). With 
equation (Bl8) 

Cl oe (C7) 


Expanding equation (Bl8) for cp = it/8, n/ 4 , and 3 «/ 8 , or obtaining 
Gi, G 2 , and G 3 in series of a’s, the sum of the G's gives 

a 2 = 2 (°-7071Gx + G 2 + O.707IG3) + a 14 - a 18 + a 30 - a3 4 (c8)^ 

The higher harmonic coefficients can be put as factors of the G n . The 
rolling-moment coefficient becomes 


c l = A 


jo.7Q71it 

J ( a l4 a 18 + a 30 “ a 34^ 

l 16 . 

L 1.2071G1 


Gi + 


it 

16 


1 + 


(a 14 a 18 + a 30 — a^) 


0 . 7071 * 

16 


1 + 


. I.207IG2 

( a i 4 — a 18 + a 3Q — a 34) 

1 . 2071 G - 




} 


(C9) 


When h n is defined as the coefficients of G n 

Ci — A (hiGi + hgG2 + h 3 Ga) 


(CIO) 


The ratio of (a 14 — a^g + a 33 — a3 4 ) to Gn can be evaluated by 
the zero-aspect-ratio theory. It is expected this ratio will not vary 
appreciably with aspect ratio. From reference 6, the loading in series 
expansion gives for equation (Bl8) 
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(T) 

\ / outboard. 


8 A 


inboard 




(cos Q sin & — M-i sin 6 cos hi 0) 


j£i 

-(-1) 2 
Hl 2 -1 


>(cn) 



-(¥) 

= 1 ' outboard 


J 


These high harmonic coefficients are small, but are not negligible for 
loading due to ailerons. The h n are tabulated for each of the cases. 


TABLE C5 


h n 

Case I 

II 

III 

IV 

<1 

VI 

hi 

0.1398 

0.1388 

0.1379 

0.1462 

0.1407 

0 . 1402 

ha 


.1963 

.1955 


.1973 

.1975 

h 3 

.1446 

.1388 

.1382 

.l4oo 

.1394 

.1397 


Spanwise Loading Distribution 


The spanwise loading distributions due to the twist distributions 
of table C4 are found at three span stations, and, since these loadings 
are not completely defined by a few terms of the assumed loading series, 
the values of loading at other span stations cannot be found accurately 
by direct use of equation (B23) and table Bl. For zero— aspect— ratio 
wings, the spanwise loading distribution due to aileron deflection is 
given at all span stations by equation (C3). The loading distribution 
for other than zero-aspect-ratio wings will fluctuate about the value 
given by equation (C3) in a manner similar to the manner that loading 
due to rolling varies about the function sin 2 <p of zero-aspect-ratio 
theory. Since the interpolation table of equation (B23) applies only to 
loadings that vary about the function sin 2 cp the loading due to 
aileron deflection can be divided by the ratio of equation (C3) to sin 2 <p 
and the resulting loading will be approximately given by sin 2<p. 
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The zero-aspect— ratio values of equation (C 3 ) can be tabulated as 

ratios of QL 5 H . . Define 
sin 2 cp 

(Cl 


E n = 


G(<Pn)/S 
sin 2% 


The .zero-aspect— ratio values of Rn can be tabulated for each aileron 
span case considered. 

TABLE C6.- Rn 



Outboard 

1 

11 

III 

0.1685 

0.4444 

0.8049 

0.1607 

.0500 

.0269 

0.2714 

.2800 

.1445 

0.3275 
• 3851 
•5119 



0.5556 0.8315 1.0000 


0.0642 0.1749 0.3358 

.1164 .3464 .3964 

.4132 .5309 .5578 


The interpolation series of equation (B 23 ) becomes 


L euk Vr l 

n=l u < 


n" _ G 


where G = - and e^ are given by table Bl. With R k tabulated 
values of loading at span stations = cos are obtained. 


TABLE C 7 .- Rv 



0.981 1/2 « 0.1738 0.2663 0.3162 0.0559 0.1484 0.3222 

.831 3/2 .1056' .2787 .3499 .0802 .2533 .3589 

.556 5/2 .0341 .2250 .4365 .2424 .4225 .4566 

.195 7/2 .0233 .1212 .5304 .6363 .7343 .7575 
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TABLE I.— ANTISYMMETRIC INFLUENCE COEFFICIENTS, Pvn> 
BEYOND TEE SCOPE OF FIGURE 1 




15.23 

16.91 

16.96 

19.25 

18.68 

21.60 

20.46 

23.91 

22.27 

26.22 

2 b . 08 

28.53 

25.86 

30.84 

j 







.NACA^ 


































































































































NACA TN 2140 


55 



<o 

>< 

Ci 

CO 

co 

«Nj 




<0 

<3 

»i 

'O" 


h 



Figure /.- Influence coefficients, p vn , for antisymmetric spanwise loading plotted as a 
function of the wing geometric parameter, H y , for values of the compressible sweep 
parameter. A# degrees. 
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Figure /- Continued. 
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Figure /.- Continued. 
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Figure /.- Continued. 
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Figure /.- Continued. 
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Figure /.- Continued . 



66 


NACA TN 2140 



NACA TN 2140 


67 



Figure /.- Continued. 
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(h) r* 3, n m 2, df 0.38/ . 
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(i) v m J, n • 3, dfO.381. 
Figure /- Concluded. 
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Figure 2.- Variation of the geometric parameter 
straight-tapered wings. 
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with taper ratio A for 
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Figure 3- Variation of lift -effectiveness parameter with aileron chord ratio, 
t • —jr- . Average trailing- edge angle about 10 * ; M *0.2. Curves from 
reference 7 . 
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(a) v =0.3827 . 


c,c 

Figure 4 - Variation of loading due-to-rolling coefficient with com- 

pressible sweep parameter degrees, for straight- tapered wings. 
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(b) v =0707! . 


Figure 4. - Continued. 
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(c) 7} =0.9239 . 


Figured - Concluded. 
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(a) \=0. 


^ c 2 p 

Figure 5 -Variation of damping-in-roll parameter with compressible 

sweep parameter degrees , for straight -tapered wings. 
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(b) A =0.25. 


Figure 5. - Continued. 
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(c) A = 0.50. 


Figure 5- Continued. 
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(d) A = 10. 


Figure 5 - Continued. 
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(e) A = 15 . 


Figure 5rConc!uded. 
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(a) ,1 - O. 


$Cj 

Figure 6- Aileron rolling-moment parameter —~,per radian, and rolling 

pc 2 " 

moment due to sideslip with dihedral — , per radian squared, for 

te av i * 

extent of unit antisymmetric angle of attack from the wing root 
outboard. 
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(b) A * 0 . 5 . 


Figure 6 - Continued 
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(c) A * i.O . 


Figure 6 - Concluded 
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Figure 7 - Rolling moment due to sideslip of wing with dihedral ” 0 

K av r 

per radian squared, for unit constant spanwise dihedral angle for 
straight - tapered wings . 
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Figure 8- Spanwise loading due to rolling of wings with various taper 
ratio and sweep angle parameters of aspect ratio parameter = 4 . 

^QV 

The curve for A-0 serves as a basis for comparison and the factor 
V 

— cp — gives the curves constant area equal to 3.454 . 

%p(A-0) 



cp(A*0) 
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Figure 8 - Concluded . 
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Figure 9 - Variation of the damping- in- roll parameter 


ratio parameter 


for various sweep angles and taper ratios. 
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0 % 




wings with various aspect ratio and sweep angle parameters. 



Figure II- Antisymmetric span wise loading distribution due to 
deflection of outboard ailerons, for c a =c, on wings hav- 
ing a taper ratio of 0.5. 





Figure 12 - Aileron rolling-moment-coefficient parameter — -r- 

'w 

radian, for outboard ailerons as a function of aspect ratio 
parameter, taper ratio, and compressible sweep parameter . 
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(c) Effect of sweep, rj g - 0.5 . 


Figure !2- Concluded . 




Figure 13.- Aileron rolling— moment parameter — — * , per radian, for — — - 4, A =0.5, as 

*av Kav 

a function of outboard and inboard aileron spans 



Figure 14 r Comparison of damping- in- roll coefficient Cj p of the 
present theory with those of the theories of references i \ 6 t and 8. 
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-JO -.15 -20 -25. -.30 -35 -40 

Experimental C lp 

Figure 15.- Correlation of theoretical and low- speed experimental damping- 
in- roll coefficient C lp for various plan forms. 










Theoretical C, 


NACA TN 2lk0 


95 



Figure !7.~ Correlation of theoretical and low-speed experimental aileron- 
effectiveness C x , per radian, due to two antisymmetrically deflected 
ailerons, for various plan forms. 
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